Convergence of Particle-Hole Expansions for the Description of Nuclear Correlations 
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The convergence properties of a multiparticle-multihole {mp-mh) configuration mixing approach 
whose purpose is to describe ground state correlations in nuclei without particle number and Pauli 
violations is investigated in the case of an exactly solvable pairing hamiltonian. Two different 
truncation schemes are tested by looking at quantities as correlation energies and single-particle 
occupation probabilities. Results show that pairing correlations present in usual superfluid nuclei 
can be accurately described using up to 6 particle-6 hole excitations, a convergence fast enough for 
envisaging extensions to fully microscopic calculations. 



I. INTRODUCTION 

One of the most powerful method of predicting nuclear 
ground state properties and excitations all over the nu- 
clear chart is the microscopic approach based on the self- 
consistent mean- field theory 1]. Within this kind of ap- 
proach, the nuclear mean- field is determined from a vari- 
ational technique of the Hartree-Fock (HF) type applied 
to a many-body Hamiltonian. In order to be tractable 
up to the heaviest nuclei, the method is usually applied 
by employing phenomenological effective two-body inter- 
actions such as the Skyrme 2] or the Gogny S] forces. 

It is well known that an accurate description of nuclear 
structure almost always requires taking into account cor- 
relations beyond the simple HF approximation. One im- 
portant class of these correlations, namely pairing corre- 
lations, are commonly treated using the BCS or Hartree- 
Fock-Bogoliubov (HFB) theories. A second class of cor- 
relations comes from the coupling of the nuclcon motion 
to collective oscillations of the mean field. Depending on 
the amplitude of collective oscillations, these correlations 
are usually derived either from the Random Phase Ap- 
proximation (RPA) theory or from the Generator Coordi- 
nate Method (GCM). Let us emphasize that, in addition 
to numerous standard nuclear properties, correlations are 
expected to play a key role in the interpretation of the un- 
usual structure recently found in some nuclei such as the 
so-called parity inversions observed in the ground state 
and low- lying excitations of exotic light nuclei 0, Q • 

Most techniques used for including correlations suffer 
from defects that may be an obstacle to an accurate de- 
scription of nuclear structure observables. For instance, 
proton and neutron numbers are not conserved in the 
BCS and HFB schemes. This deficiency is known to ren- 
der these theories inadequate for describing superfluid to 
normal phase transitions and, more generally, situations 
where pairing correlations are small. As a consequence, 
the precision of the approximations consisting in treating 
simple excitations as multi-quasiparticle states in even- 
even and odd-even nuclei becomes questionable. Simi- 
larly, the quasi-boson approximation employed in RPA 
generates violations of the Pauli principle that require 



to introduce corrections to the mean-value of one-body 
operators and render delicate a microscopic treatment of 
particle-vibration coupling. 

These difficulties can be obviated by having recourse to 
particle number projection techniques in the case of pair- 
ing or to extensions of RPA that do not make use of the 
quasi-boson approximation. However, these techniques 
complicate numerical algorithms considerably. For in- 
stance, in order to get a realistic description of weak pair- 
ing situations necessitates to perform the particle num- 
ber projections before variation [y|. In the same way, 
extensions of the RPA avoiding the violation of Pauli 
correlations as e.g., the Self-Consistent RPA Jj require 
significant additional numerical effort. 

In view of this, an alternative approach has been pro- 
posed in Ref. 8] , which consists of taking a ground state 
trial wave function in the form of a linear combination of 
multiparticle-multihole {mp-mh) operators acting upon 
the HF state. The relevance of taking as a starting point 
a set of HF single-particle states has been pointed out 
in Ref. 13 • The mp-m/i expansion is thought to be trun- 
cated at a given order and the mixing coefficients are 
determined by minimizing the total energy of the sys- 
tem. The interest of such a trial wave function is of 
course to exactly conserve nucleon numbers and be fully 
consistent with the Pauli principle. In lowest order, in 
which only Ip-lh configurations are considered, one gets 
the well-known Tamm-Dancoff approximation (TDA) |j] . 
Adding 2p-2h and higher order configurations allows one 
to describe ground state correlations of the most general 
form, including pairing and RPA correlations. In the 
simplest application of this method, a HF calculation is 
performed and the secular equation involving the matrix 
elements of the residual interaction is diagonalized. The 
coefficients of the mp-mh superposition for the ground 
state are then taken from the eigenvector belonging to 
the lowest eigenvalue. The eigenvectors belonging to the 
next few higher eigenvalues can be interpreted as approx- 
imately representing the lowest excitations of the system. 
Let us point out that the advantage of building mp-mh 
configurations from HF single-particle states lies in the 
fact that an important part of two-body correlations is 
already included in the nuclear mean field. One there- 



fore expects that the effect of the residual interaction is 
small enough to allow a fast convergence of the mp-m/i 
expansion so that only the lowest orders will have to be 
retained. 

A more elaborate version of this method is to allow 
the HF single particle themselves to be included as vari- 
ational parameters in the energy minimization. This pro- 
cedure yields additional equations from which renormal- 
ized single particle states can be derived. The method 
then is similar to the Multi-Configuration Hartree-Fock 
(MCHF) approach which is widely used in atomic and 
molecular physics [T^ • The interest of introducing such 
a consistency between the single-particle structure and 
the correlated ground state is that additional parts of 
two-body correlations are further included in the nuclear 
mean field. In this situation even lower orders in the mp- 
m/i expansion should be necessary to describe ground 
state correlations compared to the case without single- 
particle state renormalization. An approximate but nu- 
merically far less demanding way of implementing such 
a self-consistent definition of single-particle states is to 
derive them from the mean-field calculated with the cor- 
related ground state one-body density matrix. This latter 
kind of self-consistency has been tested in Ref . |3| for cor- 
relations generated by a zero-range residual interaction in 
A~180 nuclei. The result was that the Tap-mh expansion 
could be truncated at 4 particles~4 holes. 

Taking the single-particle states entering the mp-mh 
expansion from a self-consistent mean-field is the essen- 
tial difference between the method of Ref. Q and the par- 
ticle number conserving treatments of the pairing hamil- 
tonian proposed in Ref. |lll . Iia | . It differs in the same 
respect from the large-scale shell model approach which 
is widely applied to light nuclei and to heavier nuclei near 
closed shells TS]. In addition, contrary to shell model cal- 
culations, the whole single-particle spectrum is a priori 
considered in the definition of mp-mh configurations, the 
space truncation resulting from imposing a limit on the 
order a of the ap-ah configurations included. 

In principle, given a many-body hamiltonian capable 
of describing both the nuclear mean-field and the main 
kinds of correlations beyond the mean field, such an ap- 
proach can be implemented in a fully microscopic way. 
This can be done in particular by choosing a many-body 
effective hamiltonian built with the standard density de- 
pendent effective forces which have proved to give a pre- 
cise account of pairing and RPA correlations. In this 
case, the coefficients of the mp-m.h expansion can be de- 
rived by applying a variational principle to the total nu- 
clear energy. The application of the method of Rcf.Q 
then appears as a straightforward generalization of HFB 
and RPA calculations, without the above mentioned par- 
ticle number and Pauli violations. 

Still, conceivable applications to realistic situations in 
nuclei require that one is able to describe relevant cor- 
relations using a small enough number of particle-hole 
excitations in the mp-mh expansion. In order to test the 
capacity of the method to describe the correlations gen- 



erated by a pairing residual interaction, we have applied 
it to the exactly solvable hamiltonian proposed long time 
ago by Richardson [lj| . This simple model consists of a 
system of 2N identical fermions distributed over a set of 
2N equally spaced, twofold degenerate levels il5ij. The 
mp-mh method has been applied to this hamiltonian un- 
der the form of multiple pp-hh excitations, where p and 
h are the time- reversed of p and h. Particle numbers and 
pairing strengths have been explored within ranges of 
values corresponding to typical situations found in realis- 
tic nuclei. Ground state correlation energies, occupation 
probabilities and odd-even mass differences have been 
calculated and compared with the exact results given by 
Richardson for different numbers of pp-hh excitations. 
Let us emphasize that the purpose of the present work is 
not to discuss the interpretation of the calculated quan- 
tities and their connection to pairing in real nuclei. Our 
goal is simply to check the ability of the mp-mh approach 
to describe pairing correlations within a model where ex- 
act solutions exist. 

The paper is organized as follows. In Section II we 
present additional detail concerning the mp-mh configu- 
ration mixing method. We also give the Richardson ex- 
pression of the exact ground state solution of the pairing 
hamiltonian. Results obtained with the mp-m.h method 
are presented and discussed in Section III. A summary 
and conclusions are given in Section IV. 



II. THE FORMALISM 

In the mp-mh configuration mixing method the ground 
state wave function of the system is written as a super- 
position of a HF Slater determinant |0o >= TT'^/I^|0 > 

h 

and particle-hole excitations built upon it: 

I* >= Aapoh\(t>a > +'^ A\p^f^a+ah,\(t>a > + 



The creation operators associated with unoccupied and 
occupied single particle states in |0o > are denoted by 
a+ and d^ , respectively and the Aapah are configuration 
mixing coefficients to be determined. In compact form, 
eq.(^) can be written: 



M 



I* >= 5]II^U"'il'^"Pi"''i >' 



(2) 



a=0 i=l 



where i denotes an ordered set of indices specifying a 
given ap-ah configuration, \4>apiahi > is the wave func- 
tion obtained by acting with the corresponding ap-ah 
operator on |0o >= I^OpOh >i "^q is the number of con- 
figurations of ap-ah type {mo = 1) and M is an integer 
parameter. The configuration mixing coefficients A^^nah 



are obtained by minimizing the total energy of the sys- 
tem: 



d < *|i/ - A|* > _ a < ^\H - A|* > 



r)4' * 



dA 



= 0, 



(3) 



cxpah 



where A is the Lagrange parameter related to the conser- 
vation of the norm of |\1/ >. 

When H is independent of the nuclear density, the con- 
ditions Q are equivalent to the secular equation express- 
ing the diagonalization of H in the mp-mh space : 



EE 



< (Pa 



P!°'h: 



\H\ 



^"p./"'h., ^^a'pa'h — ^ ^apah (4) 



and the complex conjugate one. 

In this work we choose for H the pairing hamiltonian: 

2N 2N 2N 

H = ^ef{a^af + a+jaj) -5E E a'^ajapaf', (5) 

/ = ! /=!/' = ! 

where / denotes the time-reversed state of /. The exact 
ground state of this hamiltonian for the case of a system 
of 2N fermions has been derived by Richardson [T^ : 



N 



^^exact ~^^Y[B+\0>, 



(6) 



where the operator B^ creates a collective pair : 

N 



^t-Y.^,7^^ 



J=l 






(7) 



and [0 > is the fermion vacuum. The quantities Ei are 
the solutions of a set of N nonlinear equations [I^ : 

iV ^1 

j(#4) = l ■< 3 = 1 ■> 

and the exact ground state energy of the system is given 
by: 



N 



E = Y^E, 



(9) 



The solutions of the pairing hamiltonian can be classi- 
fied by the eigenvalues of the seniority operator: 



f 



(10) 



which counts the number of unpaired particles present in 
the system. In the mp-m/i scheme the ground state of the 
pairing hamiltonian which corresponds to seniority zero 
can be generated by considering in ^ only configurations 
ap-ah — Ppp-f3hh with even values a=2p. If one includes 
in expansion (j^J all ap-ah configurations up to M=2N, 



the wave function |^ > becomes equivalent to the exact 
wave function |vl/'=^"ct ^_ 

It is interesting to discuss the relationship between 
the mp-iah configuration mixing method presented above 
and the projected-BCS (PBCS) approach 0, [13, em- 
ployed for restoring particle number conservation in finite 
Fermi systems. In the PBCS approach the ground state 
has a pair condensate structure obtained by replacing in 
eq.(|HJ the N pair operators B^ by a unique collective 
pair : 



N 

E^^44 



(11) 



where the amplitudes Xj are determined variationally. 
The PBCS wave function can be written also in terms 
of particle- hole operators acting on the HP state [T^ : 



\PBCS >= const. 



E 

n=0 



{niy 



\HF>, (12) 



where the pair operators P/i (Pp) are formed by restrict- 
ing the summation in Hll|) to hole (particle) states and 
replacing Xj by Xp (1/xh). From (|12|l it can be seen that 
the PBCS wave function belongs to a particular subset 
of the Top-Tnh wave functions in which the mixing 
coefficients have a separable form in the indices associ- 
ated with particles and holes. It is known that, when 
the Xj's are determined from projection before variation, 
PBCS wave-functions give a satisfactory description of 
weak pairing regimes and of the crossover from normal 
to superfluid phases in Fermi systems. One therefore ex- 
pects that the more general mp-m/i form (Q is able to 
provide an accurate description of pairing correlations in 
these situations. In addition, since the Xp and l/xh are 
small, only the first few terms in (^ should be necessary. 
In order to test the convergence of expansion ^ in all 
pairing regimes, the ground state solution of the hamilto- 
nian © has been calculated using eq.|@J for two systems, 
one with 2N=8 particles distributed among 8 equidis- 
tant levels, and the other one with 2N=16 particles dis- 
tributed among 16 equidistant levels. The levels are 
twofold degenerate. These two systems simulate the level 
densities and pairing diffusivities found in weakly and 
moderately superfiuid nuclei, respectively. The constant 
g has been varied in a large range of values around typical 
pairing strengths in nuclei. 



III. NUMERICAL RESULTS 

The number of seniority-zero states of the 2N particle 
system is equal to the binomial coeflficient C^. Thus, 
for the systems with 2N=8 and 2N=16 particles the total 
number of seniority-zero states which should eventually 
be used for exactly calculating the ground states of the 
two systems is equal to 70 and 12870, respectively. We 



2N 


2p2h 


4p4h 


6p6h 


8p8h 
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16 


36 


16 
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16 


64 


784 


3136 


4900 



TABLE I: Number of app-ahh configurations for given values 
of 2a between 2 and 8 for tiie systems with 2N=8 and 2N=16 
particles. 
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64 


776 


1946 


464 












70 


64 


784 


3118 


3710 












128 


64 


784 


3136 


4862 



TABLE IL Number of configurations corresponding to various 
cut-off energies Ecut for a system of 2N particles. Left: 2N=8. 
Right: 2N=16. Ecut is in unit of the level spacing d. 



study here the convergence properties of the mp-m/i ex- 
pansion 1^ as a function of the maximum particle-hole 
order M and of an energy cut-off in the energy of the 
vap-Tnh configurations. The numbers of app-ahh config- 
urations for 2a between 2 and 8 are shown in Table ^ 
for the two systems. The highest number of configura- 
tions is reached when a is equal to half the number of 
particles. The numbers of configurations for various cut- 
off energies Ecut are given in Table E] Here and in the 
following, energies are in units of the equidistant level 
spacing d. 

The convergence of the mp-ni/i expansion towards the 
exact solution is studied below by looking at correlation 
energies, odd-even mass differences and level occupation 
probabilities. The same quantities calculated with the 
plain BCS (unprojected) approach are also given. 



A. Ground state correlation energy 

The correlation energy Ecorr in the ground state of the 
system is taken as: 



Ecorr = E{g ^ 0) - E{g = 0), 



(13) 



where E{g) is the total energy of the interacting system. 
The amount of pairing correlations in the 2N particle 
ground state can be estimated from the odd-even mass 
difference defined by Q : 

P{2N) = 2E{2N - 1) - E{2N) - E{2N - 2), (14) 

where E{N) is the ground state energy for N fermions. 
In the calculations we adopted the convention of Ref . [l5l| . 



i.e., for a system of TV particles {N odd or even), the 
pairing force is effective among N doubly-degenerated 
levels. Notice that P{2N) = for non-interacting par- 
ticles {g — 0). Therefore, the mean- field effect which is 
contained in odd-even mass differences for real nuclei [ig 
is not present in the equidistant model used here. Thus, 
P{2N) can be taken as a genuine measure of the intensity 
of pairing correlations. 

Assimilating P{2N) with the pairing gap A, the 
weak (strong) pairing regime corresponds to values much 
smaller (larger) than 1 of the ratio rj — A/d. For de- 
formed nuclei in the rare earth region A is about 1 MeV 
while d is of the order of 400-500 keV. Thus, for these 
nuclei rj is of the order of 2-3. As we shall see below, 
for values of 77 in this range the correlation energies eval- 
uated in the BCS approximation are largely underesti- 
mated while the mp-m/i expansion converges quickly to 
the exact results. Let us add that, in units of d, the con- 
stant g corresponding to typical nuclear pairing strengths 
is about 0.5. 
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FIG. 1: Percentage error on correlation energies as a function 
of the pairing strength g for the 2N=8 and 2N=16 cases. The 
curves correspond to the results obtained by cutting the mp- 
mft expansion at 2p2h, 4p4h and 6p6h. 



We first analyze the convergence of the mp-m/i expan- 
sion for the two quantities Ecorr and P{2N) when in- 
cluding successively all 2p2h, ApAh, . . . , configurations in 
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0.2 
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-0.924 


0.36 
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-0.927 


0.00 


0.303 
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0.00 
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0.00 
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0.00 




8 


-4.136 
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0.00 
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2.81 




6 
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0.05 
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0.00 
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A£<:orr(%) 


P{2N) 


AP(2iV)(%) 


0.18 
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-1.659 


0.64 


0.305 
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-1.669 


0.00 


0.318 


0.00 




8 


-1.669 


0.00 


0.318 


0.00 
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0.00 


100.00 


0.00 


100.00 




Exact 


-1.669 




0.318 
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-6.883 


20.92 


1.941 


49.51 




4 


-8.250 


5.22 


3.301 
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-8.652 
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3.781 


1.64 
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-8.702 


0.03 


3.842 


0.05 
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-6.636 


23.76 


2.580 


32.88 




Exact 


-8.704 




3.844 
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-9.149 


29.37 


2.598 


57.16 
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-11.711 


9.59 


4.746 


21.74 
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-12.738 


1.66 


5.819 


4.04 
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-12.938 


0.12 


6.047 


0.28 




BCS 


-10.357 


20.05 


3.701 


43.96 




Exact 


-12.954 




6.064 





TABLE III: Correlation energies Ecorr and odd-even mass 
differences P{2N) for a system with 2N=8 particles. The 
results are shown for several values of the pairing strength g 
and with the mp-mh expansion cut at various orders. The 
BCS and exact results are also shown. In the BCS case, the 
value given in the column P(2N) is the one of the pairing gap. 



the ground state wave function (2) (without any cut-off). 
The resuhs of the calculations for the systems with 8 and 
16 particles are shown in Tables ITlIlI VI and in Fig^ The 
errors for correlation energies, AEcorr, and for the odd- 
even mass differences, AP, are calculated relative to the 
exact values. As mentioned above, for the system with 
2N=8 particles an expansion up to 8p — 8h corresponds 
to the exact solution. 

In Tables lnillVl the results are shown for several values 
of the pairing strength g. In the BCS approximation 
there is a sharp transition between the normal and the 
superfluid phase which appears at gc ~ 0.31 for 2N=8 
and gc ~ 0.24 for 2N=16. Below these critical values 
of the pairing strength the correlation energies are zero 
in the BCS approximation. This is not the case for the 
exact solution, which shows a rather smooth increase of 
the correlation energy with the pairing strength. 

From Tables IIIIIIVI we also notice that, in the BCS 
approximation, the correlation energies are strongly un- 
derestimated. The error made is more than 20% for all 
the values of the pairing strength listed in Tables llllt 
IIVI The same large differences are obtained between the 
pairing gaps associated in BCS with the odd-even mass 
differences, and the exact values of P(2N). It is worth 
stressing that the errors in the BCS correlation energies 
remain rather large, more than 10%, even for large values 
of the pairing strength for which the ratio 77 between the 



TABLE IV: Same as Table Ell for 2N=16 particles. 



BCS pairing gap and the level spacing (which in our case 
is the energy unit) is of the order of 20. 

Let us emphasize however that in BCS or HFB calcu- 
lations applied to real nuclei, the pairing strengths for 
neutrons and protons - or the matrix elements of the 
pairing residual interaction - are adjusted in order to 
correctly describe pairing correlations in strongly super- 
fluid nuclei. Consequently, BCS or HFB results should 
be taken with caution only in weak pairing regimes. 

From Tables llilll VI one can see that the mp-mh expan- 
sion is converging rapidly to the exact results. Thus, for 
the system with 16 particles a truncation to 6p-6h pro- 
vides rather accurate values for the correlation energies 
{AEcorr ^ !%)• This truncation order should be com- 
pared to the maximum possible order of the mp-mh ex- 
pansion, which is 16p-16h in this case. As seen in TableQ] 
for a truncation at 6p-6h the number of configurations is 
equal to 3136, while the total number of seniority-zero 
configurations is 12870. 

Next, we examine the convergence of the results as 
a function of the energy cut-off Ecut ■ In this truncation 
scheme we consider only mp-mh configurations whose ex- 
citation energies are smaller than Ecut- The number of 
configurations for various Ecut values are shown in Ta- 
ble |n] In Fig[21 the evolution of AEcorr with Ecut is 
shown for several values of the pairing strength g. 

From these figures and from Table HTl one can see that, 
for the values of g of physical interest {g r^ 0.5) the num- 
ber of configurations needed to achieve the same accu- 
racy as previously is smaller. Thus, in the case of the 
system with 16 particles and for g=0.54 one needs a 
cut-off energy Ecut=5Q in order to get an accuracy of 
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FIG. 2: Percentage errors on correlation energies as a function 
of cut-off energy Ecut- The curves corresponds to various 
pairing strength values listed in the inset. The results are for 
the systems with 2N=8 and 2N=16 particles. 



AEcorr ~ 0.5%. As seen in Table [^ for this cut-off en- 
ergy one selects 3250 configurations up to 8p-8h. On the 
other hand, in the truncation scheme based on the or- 
der of the particle-hole expansion, to achieve the same 
accuracy one needed to consider all configurations up to 
6p-6h that is, from Table HJ 3984 configurations. This is 
about 20% more than in the truncation based on cut-off 
energy. In particular a large number - about one third 
- of the 6p-6h configurations of Table Q] are eliminated 
by the cut-off. Since the same precision is achieved, this 
means that these eliminated configurations do not con- 
tribute to the correlation energy. Let us mention that a 
50 MeV cut-off in the excitation energy of scattered pairs 
is consistent with the results of BCS or HFB calculations 
in actual nuclei, even in strongly superfluid ones. 
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FIG. 3: Occupation probabilities v„ for the single- 
particle levels n. The calculated points are joined 
by straight lines. The curves correspond to different 
truncation orders specified in the insets, the full line 
showing the exact values. The results are for the sys- 
tem with 2N=8 particles. 



B. Occupation probabilities 

The occupation probability v'^ of the single-particle 
level n is given by : 



vl^^<'i'\atan + a+an\'i'> . (15) 

The occupation probabilities of all single-particle levels 
in the cases 2N=8 and 2N=16 are plotted in FigslSHfor 
three values of the pairing strength g. It can be seen 
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FIG. 4: Same as Figglfor 2N=16 particles. 



that for small values of g the exact values of occupation 
probabilities are rather well reproduced by a truncation 
to 2p2h configurations. 

However, for the highest values of <? listed in Tables UTTl 



IVli.e., g=l.O for 2N=8 and g^OM for 2N=16, one needs 
to introduce up to 6p-6h configurations in order to get 
an accurate description of occupation probabilities in the 
mp-mh scheme. This truncation order appears consistent 
with the behaviour of correlation energies discussed in the 
previous subsection. 



IV. SUMMARY AND CONCLUSIONS 

In this paper we have analyzed the convergence prop- 
erties of the multiparticle-multihole {mp-mh) configura- 
tion mixing approach. The study was done for a pairing 
hamiltonian with a constant pairing force, which can be 
solved exactly. For the single-particle model we have cho- 
sen a sequence of doubly-degenerate and equidistant lev- 
els, half-filled with one kind of fermions. The parameters 
of the model and the number of particles have been cho- 
sen so as to simulate the physical situation met in light 
and medium nuclei. The pairing strength g has been var- 
ied from zero up to about twice the value corresponding 
to typical conditions in atomic nuclei. 

We have shown that a truncation based on the mp-mh 
expansion is converging rather rapidly to the exact results 
for the ground state correlation energies and occupation 
probabilities of single-particle levels. For instance, with 
a value of g corresponding to standard nuclear pairing, 
one can get an accuracy of about 1% for the correlation 
energy in a system of 16 particles if one considers config- 
urations up to 6p-6h. The number of these configurations 
is about four times smaller than the total number of pos- 
sible mp-mh configurations with seniority zero that the 
system can form. An even faster convergence without 
precision loss is found by selecting configurations whose 
excitation energies are limited by a cut-off correspond- 
ing to the maximum energy of scattered pairs. These en- 
couraging results indicate that applications of the mp-mh 
approach to fully microscopic self-consistent calculations 
in nuclei should be possible with present-day numerical 
capabilities. 
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